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Abstract 

Design and analysis of the inlet for a 
rocket based combined cycle engine is discussed. 
Computational fluid dynamics was used in both the 
design and subsequent analysis. Reynolds averaged 
Navier-Stokes simulations were performed using both 
perfect gas and real gas assumptions. An inlet design 
that operates over the required Mach number range 
from 0 to 1 2 was produced. Performance data for 
cycle analysis was post processed using a stream thrust 
averaging technique. A detailed performance database 
for cycle analysis is presented. The effect of vehicle 
forebody compression on air capture is also examined. 

Nomenclature 

A cross sectional area 

D drag 

F stream thrust 

H total enthalpy 

M Mach number 

Q dynamic pressure 

R gas constant 

U one-dimensional velocity 

V velocity vector 

c specific heat at constant pressure 

f tlow' direction unit vector 

h enthalpy 

th massflow 

p pressure 

r radius 

x distance from station 3 

p density 
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subscripts 

0 station 0, freestream 

1 station 1 

2 station 2 

3 station 3 

b body 

c capture 

cl cowl 

cb centerbody 

d diverter 

t total condition 

superscripts 

stream thrust averaged quantity 


Introduction 

One of the three primary goals of NASA's 
Aeronautics and Space Transportation Technology 
program is to enable low cost access to space by 
developing advanced space transportation concepts 
and technologies. The key to reducing space launch 
costs is developing a reusable vehicle with a short turn- 
around time. Highly reusable implies a very robust 
Single-Stage-to-Orbit (SSTO) vehicle. NASA's 
Glenn Research Center has undertaken a program to 
demonstrate such a vehicle. The concept, called the 
Trailblazer 1 (figure 1), is powered by a Rocket Based 
Combined Cycle (RBCC) engine (figure 2). This 
engine is designed to operate efficiently from lift-off 
to orbit by integrating a rocket and ramjer. The 
system combines the high thrust-to- weight/low speed ic 
impulse characteristics of the rocket with the high 
specific impulse/low thrust to weight characteristics ot 
the ramjet. The engine operates in air-breathing modes 
from lift-off to between Mach 10 and 12, at which point 
the airbreathing engine flow path is closed ott, and the 
rocket is turned back on to carry the vehicle out ot the 
atmosphere and into orbit. Three semi-circular engine 
pods are located near the aft portion of the vehicle. 

The pods allow for diversion of the boundary layer, 
simplify centerbody actuation and sealing, and enable 
integration of the nozzles with the vehicle base. 
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An important component of the RBCC engine 
is the inlet. This device must efficiently capture and 
compress the air over the entire range of air-breathing 
Mach numbers. The objective of this work was to 
design an inlet and provide detailed performance 
data for engine cycle analysis. Computational fluid 
dynamics (CFD) was used in both the design and 
analysis process. CFD provided a fast and efficient 
method for obtaining high fidelity performance data 
for the Trail blazer design effort. 

Inlet Design 

Traditional design methods involve point 
designs at a cruise condition and are not applicable 
to this problem. The resulting inlets have very good 
cruise performance but may suffer at off-design 
conditions. Inlets for airbreathing SSTO vehicles must 
trade high efficiency at one Mach number for good 
performance over the entire Mach range. Weight and 
wetted area must also be considered. Several design 
decisions were initially made to guide the process. A 
mixed compression inlet with a translating centerbody 
was chosen. Moving the centerbody allows the inlet 
to start and varies the contraction ratio over the Mach 
number range. For rocket only operation (in vacuum 
conditions), the inlet duct must be completely closed 
off. The inlet duct was designed to accomplish this by 
fully retracting the centerbody. The inlet is intended 
to operate without boundary layer bleed to keep the 
system weight to a minimum and maintain simplicity. 

The key locations for cycle analysis are 
indicated in figure 2. Station 0 is the freestream. 
Station 1 is located at the inlet spike tip. Station 2 is 
the inlet throat and is used as the nominal location for 
the beginning of supersonic combustion/fuel injection. 
Station 3 is at the end of a fixed hub over which the 
centerbody translates. The inlet geometry and key 
design parameters are shown in figure 3. The cowl lip 
radius is used to size the inlet. The maximum radius 
of the centerbody was chosen so that the annular area 
formed by the cowl lip and centerbody is equal to the 
area at station 3. A 1 2-degree cone angle was chosen 
to minimize length. At station 2, the inlet throat is 
angled at 15 degrees towards the axis to minimize 
length and wetted area in the diffuser/scram combustor 
portion of the flowpath. A back step in the centerbody 
is placed here to isolate the forward portion of the 
inlet from back pressure feeding forward through the 
subsonic portion of the boundary layer and can serve 
as a location for fuel injectors. A constant area section 
to facilitate scramjet ignition is located aft of the step. 
The angle at the end of the centerbody is specified at 
20 degrees. A Hat section is located on the forward 


portion of the cowl to allow the inlet to be oversped 
(shock inside the lip). 

Isentropic inlet contours for Mach 6 operation 
were generated using a method of characteristics design 
code' based on the aforementioned constraints. Mach 6 
was chosen as the key point in the flight envelope. 
Beyond this point, the required inlet contraction ratio 
does not change appreciably with Mach number. 
Centerbody shock on lip was specified at this Mach 
number and the shoulder on the centerbody was placed 
to cancel the reflected cowl shock. Beyond Mach 6, 
the shock angle does not change significantly with 
increased Mach number due to the hypersonic Mach 
number independence principle. Thus, the reflected 
shock will remain near the shoulder for all hypersonic 
Mach numbers. 

Several perfect gas CFD analyses using the 
NPARC code were done on preliminary configurations. 
The CFD solutions were used to evaluate the designs 
and provide guidance in adjusting the key geometric 
parameters. In particular, predictions of shock strength 
and location, potential boundary layer separations, and 
throat Mach number were used to adjust the inlet lines. 
The final inlet geometry is completely described in 
table I. 

CFD Method 

The flowfield was assumed to be axisymmetric 
so that two-dimensional CFD could be used. This 
assumption is valid for a large majority of the inlet 
flow. Only the effects of the inlet sidewalls/endwalls 
would alter the flowfield in the circumferential 
direction. The assumption significantly reduces the 
number of grid points required to describe the geometry 
and hence reduces the computational cost for a solution. 
Two-dimensional CFD yielded timely results for the 
inlet design process and enabled a large number of 
cases to be run for cycle analysis. Solid walls were 
specified with an adiabatic no-slip wall boundary 
condition. 

Grid Generation 

The inlet geometry is defined by a series of 
cubic splines and other simple geometric shapes. The 
spline data was used to create a series of points on the 
inlet surface that were read into the grid generator. 

The grid generation software used for this 
project was GRIDGEN 4 . It is an interactive software 
package with a user-friendly graphical interface. It 
is capable of producing high quality structured and 
unstructured grids for complex two and three- 
dimensional geometries. 
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All grids were generated such that the first 
grid point away from a viscous wall was placed at an 
average inner variable distance, y-K of 1 . This distance 
was computed based on assumed average skin friction 
and flowfield values, and was subsequently checked for 
accuracy. It has been shown that proper resolution of 
the boundary layer is not possible unless at least one 
grid point is within a y+ of two\ 

A typical grid is shown in figure 4. The grid 
uses a blocked structure with point-to-point matching 
at the block interlaces. For started supersonic cases, it 
was not necessary to model the area above the cowl lip 
because the flow in this region does not influence the 
region of interest and is not captured by the inlet. Thus, 
the grid above the cowl lip was not used for these cases. 
A grid sensitivity study was performed with the GASP 
code. Doubling the number of grid points in each 
direction did not significantly affect the solution. 

NPARC 

The NPARC code* was used for the perfect 
gas analyses. NPARC is developed and supported 
by the NPARC Alliance, a collaboration between the 
NASA Glenn Research Center and the United States 
Air Force’s Arnold Engineering Development Center. 
NPARC solves the Reynolds averaged Navier-Stokes 
equations in strong conservation law form using the 
Beam-Warming approximate factorizations algorithm 7 . 
Spatial discretization is performed using a central 
difference scheme. Jameson style artificial dissipation^ 
is added for stability and to smooth shock oscillations 
and odd-even grid point decoupling. The code uses a 
perfect gas equation of state. The code has several 
options for modeling turbulence varying from algebraic 
zero transport equations to one and two equation 
models. Both the Chien low Reynolds number k-£ y and 
Wilcox’s k-o) ln models were investigated for this work. 
It was found that the two-equation models produce very 
similar answers. All data reported from the NPARC 
code was generated using the k-e model. 

GASP 

The GASP code 11 was used for the real gas 
analyses. GASP is a commercially available code that 
is developed and sold by Aerosoft Inc. For this analysis 
the code solved the Reynolds averaged Navier-Stokes 
equations using a third-order accurate upwind scheme 
and Roe’s flux difference splitting 12 . Time 
discretization was done by a two-factor alternating 
direction implicit (ADI) scheme. Van Albada limiters 1 ^ 
were employed for stability in the presence of 
discontinuities. The chemical kinetics model ot Kang 


el. al. 14 was used. The equilibrium vibrational 
thermodynamics model was selected. The k-ou 
turbulence model was found most robust and was used 
for all GASP results herein. To ensure that there were 
no discrepancies between codes due to the difference in 
turbulence models, both NPARC with the k-e model 
and GASP with the k-co were run using the perfect gas 
assumption at several Mach numbers. Integrated results 
agreed within one percent for all Mach numbers except 
Mach 6, where the difference was approximately 5 
percent in total pressure recovery and 2 percent in 
Mach number. 

Inlet Performance Analysis 

Results were obtained for freestream Mach 
numbers from 0.5 to 12. The corresponding station 1 
Mach numbers, used for the inflow boundary condition, 
were derived based on a 10-degree conical shock and 
re-expansion to freestream pressure (table 2). A 
subsequent two-dimensional analysis of the vehicle 
forebody provided more accurate results and is 
discussed in a later section. For the subsonic cases 
(Mach 0.5 and 0.8), the centerbody was completely 
extended and the mass flux on the outflow boundary 
was regulated to control the inlet flow. For supersonic 
cases up to Mach 6 (the transition Mach number from 
subsonic to supersonic combustion), several super- 
critical (no back pressure imposed) cases were run to 
determine a centerbody position which provided the 
maximum contraction ratio. Then using this centerbody 
position, a back pressure was applied through the exit 
boundary condition to simulate the effect ol the 
combustion process. This back pressure was increased 
in several increments until the inlet was unstarted. The 
data at maximum back pressure is presented here. 
Beyond Mach 6, only super-critical cases were 
necessary due to the supersonic combustion process. 
Three inlet contraction ratios were run at these Mach 
numbers to provide a range of data (figure 5). Billig’s 
suggested contraction ratio 1 * 1 is shown for comparison. 
The maximum contraction ratio was found to be about 
16 for all hypersonic Mach numbers (Mach 6 and 
above). Results using the perfect gas code (NPARC) 
were run at all Mach numbers. Real gas results (GASP) 
were obtained for Mach 4 and higher. 

Representative Mach number contours, from 
the NPARC code, are shown in figure 6 for Mach 0.8, 
Mach 2, and Mach 6 (started) operation. A mildly 
supercritical solution is presented for Mach 0.8. 
Choking at the centerbody shoulder determines the 
maximum airflow. A normal shock is located just 
downstream. The plot also shows a region of low 
pressure on the cowl, which creates a suction force. 
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At Mach 2, a classical lambda-shock pattern appears 
just forward of the cowl lip. This case is slightly sub- 
critical, 10 percent more airflow was observed at very 
low back pressure. The shock induced separation 
reduces the How area and the flow chokes at the 
aerodynamic throat. The Mach 6 contours show shock- 
on-lip operation. It is evident that the cowl shock is not 
completely canceled on the forebody at this Mach 
number. Viscous effects were not accounted for in the 
MOC design code and alter the Mach number at which 
the shock is canceled. Because the inlet is not a point 
design, the exact Mach number where the shock is 
canceled is not significant. For all started cases the 
cowl shock wave strikes the centerbody shoulder and 
the resulting pressure gradient separates the boundary 
layer. This separation is the dominant mechanism in 
the flowfield. Centerbody surface static pressure 
distributions for Mach 6 are presented in figure 7. 
Several back pressured cases, super-critical to 
maximum back pressure (350p 0 ), are shown. 

Stream Thrust Averaging 

An averaging procedure which conserves 
mass, momentum and energy was used to compute 
integrated one-dimensional performance data from the 
two-dimensional CFD solutions. This stream thrust 
averaging (STA) technique provides a consistent set of 
data for use in cycle analysis. In STA the distorted 2-D 
profile is passed through an idealized control volume 
where the flow' is completely mixed at the exit. The 
equations are derived by applying mass, momentum, 
and energy conservation to this process. The STA 
technique for a perfect gas is presented. To simplify 
the procedure, a unit vector in the How direction, / , is 
chosen based on the geometry. Then, for each plane 
of two-dimensional data the following integrations are 
performed. 


Stream thrust. 

F= f (pVdA)Vf+\pdA-f 

A A 

Massflow, 

m= \pV il A 

A 

Stagnation enthalpy. 

If-- 

H = — (pV-dA)(h + u -) 
hi \ 2 


These three quantities are substituted into the following 
quadratic equation that is solved for the stream thrust 
averaged velocity, U . 


r F- HR „ 

U(] ) U + -0 

2 r ht c„ 


Selecting between two resulting solutions for U is done 

by applying the second law of thermodynamics. The 

remaining variables can then be found as follows. 

- hi _ F r c p 

p-^ y p- p(J\ 1 1=-^- 

y UA A y pR 

A similar technique for the chemically reacting flow 
solutions was used. 16 


Mass Capture 

Inlet mass capture versus station I Mach 
number is plotted in terms of the area ratio A/A, in 
figure 8. Subsonic, supersonic external compression, 
and started regions are denoted by different symbols. 
The subsonic and external compression cases were run 
with a low back pressure specified at the outflow 
boundary to determine the maximum airflow possible. 
Sub-critical cases were run by specifying a lower mass 
flux at the outflow boundary. The inlet starts between 
Mach 2 and 2.5 and the shock-on-lip condition occurs 
at Mach 6. The lines connect the critical (maximum) 
airflow points. 

Recovery 

Stream thrust averaged total pressure recovery 
is presented in figure 9. For the perfect gas solutions, 
isentropic perfect gas relations were applied to the 
stream thrust averaged pressure and Mach number 
to obtain the stagnation condition. The reference 
stagnation condition, at station 1 , was also obtained 
from isentropic relations. For the real gas solutions, the 
equilibrium air stagnation condition was determined 
from the stream thrust averaged state using the HAP 
code. 17 The stagnation condition of the reference state, 
station 1, was also calculated this way. Recovery was 
independent of massflow for the subsonic and 
supersonic external compression cases. Between Mach 
4 and 6 the combustion process transitions between 
subsonic and supersonic combustion modes so 
recoveries at both stations 2 and 3 are presented. Both 
perfect gas and real gas station 2 recoveries are shown. 
Mil. E-5007D is shown for comparison. The spike 
positions corresponding the recovery data are presented 
in table 2. 


Inlet Drag 

In the Trailblazer force accounting system, 
inlet drag consists of the force on the forward facing 
cowl surface and the additive drag due to spillage. 
The cowl force is computed by simply integrating 
the pressure and shear stress over the surface area. 
Additive drag is defined as the pressure force exerted 
on the captured stream tube. It is computed by 


4 

American Institute of Aeronautics and Astronautics 



AIAA 99-2239 


applying the control volume pictured in figure 10 to 
the CFD solutions. Figure 1 1 shows the drag forces as 
a function of station 1 Mach number. Lines connect the 
critical airflow points. At subsonic flight conditions the 
cowl lip pressure is less than ambient and helps offset 
the additive drag. At these Mach numbers the drag is 
the lowest at the critical mass capture. For supersonic 
cases below Mach 2.5 (external compression), 
separation and choking just aft of the cowl lip forces 
the normal shock onto the forward portion of the 
centerbody and the additive drag is large. At Mach 2.5, 
the inlet starts and the additive drag is significantly 
reduced. 

Comparison Between Real and Perfect Gas 
Solutions 

Figure 1 2 compares real and perfect gas 
solutions at station 2. The percentage difference in the 
data between perfect and real gas solutions increases 
with Mach number. For recovery, good agreement is 
obtained up to Mach 6. At Mach 10, the perfect gas 
recovery is 30 percent higher than the real gas value. 

For throat Mach number, the real gas result is just over 
10 percent higher than the perfect gas result at Mach 6 
and increases to 20 percent at Mach 10. 

Vehicle Forebody Precompression 

The Trail blazer vehicle differs from many 
other SSTO and hypersonic vehicles in the fact that 
the vehicle forebody is not designed as a compression 
surface for the inlet. The Trailblazer design has chosen 
to trade the efficiency of forebody compression for the 
simplicity and higher structural and volumetric 
efficiency of an axisymmetric vehicle. It was assumed 
that the captured airflow was processed by the shock 
resulting from a 10-degree cone, then re-expanded to 
freestream pressure at station 1. A computational study 
was performed to check the validity ol this assumption. 

An axisymmetric perfect gas CFD analysis, 
using NPARC, was performed on the forebody at Mach 
2, 4, 6, 8, 10, and 12. The grid contained 239 nodes in 
the axial direction and 245 nodes in the radial direction. 
Increasing the grid size in both directions by 50 percent 
did not significantly affect the solution. The 
axisymmetric solution yields flowfield variables as 
a function of radial distance, r, tor a given x-location. 
For each axial location, the solution was then mapped 
onto the projected aerodynamic capture area of the 
inlet. This capture area is offset from the vehicle body 
to allow for the boundary layer diverter (figure 13). 

The boundary layer diverter radius, r d , is 1 1 1 percent of 
the vehicle body radius, r h . Stream thrust averaging of 


the flowfield over this area was then performed. A plot 
of mass capture versus axial location is presented in 
figure 14. The increase in mass capture over that of a 
freestream inlet is attributed to forebody compression. 
This increase in mass capture is reduced asymptotically 
as station I is moved aft and the How re-expands to 
nearer freestream pressure. The reduction in capture 
for the high Mach numbers at forward locations on 
the body is due to the forebody shock lying inside the 
projected capture area at station 1 . When this occurs, 
a portion of the captured area is not yet compressed by 
the forebody shock. However, it is important to note 
that all the flow in the captured stream tube can still be 
compressed by the shock before entering the inlet. In 
other words, the forebody shock lies outside the cowl 
lip. The location of station 1 for the current Trailblazer 
configuration is represented by the vertical dashed line. 
Data comparing the effect of the forebody 
precompression to the original assumption of isentropic 
expansion back to freestream pressure are shown in 
figure 15. At Mach 10 the pressure is twice that of 
freestream. Mass capture is likewise affected. A 70 
percent increase in mass capture is seen at Mach 10. 

As shown in figure 15c, the stream thrust averaged 
recovery is also increased due to curvature of the 
forebody shock. Stream thrust averaged station 1 Mach 
number plotted versus freestream Mach number in 
figure I5d shows a reduction in station 1 Mach number 
from the original assumption. Concomitant with these 
effects is a radial variation in flow properties from the 
diverter surface to the bow shock. The effect of this 
distortion on inlet performance and operability is 
currently unknown. 

Conclusions 

Two-dimensional computational fluid 
dynamics was used to design and analyze the inlet tor 
the Trailblazer vehicle. The inlet is designed to operate 
from Mach 0 to 1 2 without boundary layer bleed. A 
maximum contraction ratio of 16 tor all hypersonic 
Mach numbers was obtained. Both perfect and real 
gas assumptions were used in the CFD analyses. The 
perfect gas and real gas solutions were in agreement 
up to Mach 4 for averaged throat Mach number and up 
to Mach 6 for recovery. At Mach 10, the perfect gas 
solutions yield recoveries 30 percent too high and 
Mach numbers 20 percent too low. Mass capture, total 
pressure recovery and inlet drag are presented for the 
entire Mach range. 

An analysis of the integration of the inlet with 
the vehicle forebody was also undertaken. The results 
indicate that the compression due to the vehicle 
forebody results in higher static pressure, mass capture 
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and recovery. At Mach 10 the station 1 mass capture is 
increased by 70 percent, the pressure is increased by a 
factor of 2. and the recovery is increased by 20 percent. 
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Table 2. Centerbody spike translation schedule 





Figure 1 . Trailblazer vehicle 
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Trailing Edge of 



Figure 2. Trailblazer RBCC engine 



Figure 3. Inlet geometry definition 
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b. Reduced grid for started supersonic cases 
Figure 4. Computational grid 
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Freestream Mach Number, M 0 

Figure 5. Contraction ratios analyzed 
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a. M n = 0.8 

Figure 6. Mach number contours 
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Figure 1 1 . Inlet drag 



Figure 1 2. Percent difference between real and perfect gas recoveries and Mach number at station 2 
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Cross-section 



Figure 1 3. Projected capture area tor precompression analysis 



Figure 14. Variation of mass capture with station 1 location 
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a. Pressure 



b. Mass capture 

Figure 15. Effect of forebody compression 
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Freestream Mach Number, M 0 

c. Recovery 



Freestream Mach Number, MO 

d Mach number 
Figure 15. Continued 
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